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Modeling Quasi-elastic Form Factors for Electron and Neutrino Scattering 

H. Budd a , A. Bodek a and J. Arrington b 
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We calculate the total and differential quasielastic cross sections for neutrino and antineutrino scattering on 
nucleons using up to date fits to the nucleon elastic electromagnetic form factors G V E , GJ, G P M , GJf, and weak 
and pseudoscalar form factors. We find that using the updated non-zero value of has a significant effect on 
both the total and differential neutrino and antineutrino quasielastic cross sections. Previous extractions of the 
weak axial form factor from neutrino scattering data are sensitive to the assumptions that were used for the vector 
form factors. We perform a re-analysis of previous neutrino data using updated form factors and obtain updated 
value of the axial vector mass. (Presented by Howard Budd at Nulnt02, Dec. 2002, Irvine, CA, USA pQ) 



1. INTRODUCTION 

Experimental evidence for oscillations among 
the three neutrino generations has been recently 
reported j2j. Since quasielastic (QE) scattering 
forms an important component of neutrino scat- 
tering at low energies, we have undertaken to in- 
vestigate QE neutrino scattering using the latest 
information on nucleon form factors. 

Recent experiments at SLAC and Jefferson Lab 
(JLab) have given precise measurements of the 
vector electromagnetic form factors for the proton 
and neutron. These form factors can be related 
to the form factors for QE neutrino scattering by 
conserved vector current hypothesis, CVC. These 
more recent form factors can be used to give bet- 
ter predictions for QE neutrino scattering. 

2. EQUATIONS FOR QE SCATTERING 

The hadronic current for QE neutrino scatter- 
ing is given by |3] 

<p(p 2 )|J+ |n(pi) >= 



lxF i {q * ) + m»lMMl 



2M 

2\ , q\lbFp(q 2 



M 



u(pi), 



where q = k u — k^, £ = (fi p — 1) — /i„, and M = 
(m p + m n )/2. Here, /i p and ji n are the proton 



and neutron magnetic moments. We assume that 
there are no second class currents, so the scalar 
form factor F v and the tensor form factor F\ 



need not be included, 
the cross section is 



Using the above current, 



da"' v 



M 2 G 2 F cos 2 9 c 



dq 2 
A{q 2 ) T 

where 

A(q 2 ) 
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M 



M 
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4M Z 



4q 2 ReF v *tF 2 
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B(q 2 ) - -j^ReFX(F v 
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\Fa\ 
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Although we have have not shown terms of order 
(mi/M) 2 , and terms including Fp(q 2 ) (which is 
multiplied by (mi/M) 2 ), these terms are included 
in our calculations 3 .) The form factors F v (q 2 ) 
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9a 


-1.267 


Gp 


1.1803xlCT 5 GeV~ 2 


cos 9 C 


0.9740 


M P 


2.793 ^jv 




-1.913 fj, N 


e 


3.706 /zjv 




0.71 GeV 2 



Table 1 

The most recent values of the parameters used 
in our calculations (Unless stated otherwise). 



and £,Fy(q 2 ) are given by: 



F v {q 2 ) = 



4M 2 



1 - 



4AP 
GW)~GUl 2 ) 



1 



4AT 



We use the CVC to determine G\{q 2 ) and 
G\ I (q 2 ) from the electron scattering form factors 
G E (q 2 ), G E (q 2 ), G P M (q 2 ), and G^q 2 ): 

G E (q 2 ) = G p E (q 2 )-G E (q 2 ), 

G V M {q 2 ) = G P M {q 2 )-G n M {q 2 )- 

Previously, many neutrino experiment have as- 
sumed that the vector form factors are described 
by the dipole approximation. 



G D (q 2 ) = 



1 



My 



G D (q 2 ), 



M v = 0.71 GeV 2 



G% = 0, 



u m — 



/J, P G D (q 2 ), Gm = n n G D (q 2 )- 



We refer to the above combination of form factors 
as 'Dipole Form Factors'. It is an approximation 
that is improved upon in this paper. We will refer 
to our updated form factors as 'BBA-2003 Form 
Factors' (Budd, Bodek, Arrington). Tablenjsum- 
marizes the most up to date values of the coupling 
constants and magnetic moments that we use in 
our calculations. Note that G E , G P M , and G E are 
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Figure 1. Our fits to G e /Gd, using cross section 
data only (solid) , and with both the cross section 
and polarization transfer data (dashed). The di- 
amonds are the from Rosenbluth extractions and 
the crosses are the Hall A polarization transfer 
data. Note that we fit to cross sections, rather 
than fitting directly to the extracted values of G P E 
shown here. 



positive, while G^ and the axial form factor Fa 
are negative. 

The axial form factor is given by 

9A 



F A (q 2 



1 



Mi 



This form factor needs to be extracted from QE 
neutrino scattering. However, at low Q 2 , this 
form factor can also be extracted from pion elec- 
troproduction data. 

Previous neutrino experiments used gA—— 1-23, 
while the best current value is —1.267. The 
world average from neutrino experiments for Ma 
is 1.026 ± 0.020 GeV 0]. The value of M A ex- 
tracted from neutrino experiments depends on 
both the value of gA and the values of the electro- 
magnetic form factors which are assumed in the 
extraction process. Since we are updating these 
form factors, new values of Ma are extracted from 
previous neutrino data using the better known 
values for 17^4 and the vector form factors. 

Ma can also be determined from pion electro- 
production, which yields a world average value of 
1.069 ± 0.016 GeV 0. This value should be re- 
duced by 0.055 GeV when compared to Ma as 
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Figure 2. Our fits to G p m /^, v Gd- The lines and 
symbols have the same meaning as Figure Q] 
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Figure 3. Our fit to GJj/ ^ n Go- The lines and 
symbols have the same meaning as Figure ^ 



measured in neutrino data because of additional 
corrections g]. Therefore, pion electroproduction 
experiments predict that Ma should be 1.014 ± 
0.016 GeV in neutrino scattering. 

In this communication, we show that the value 
of 1.026 ± 0.02 GeV g] as measured from the av- 
erage of all neutrino scattering should also be re- 
duced by 0.025 GeV to account for incorrect vec- 
tor form factors used in the past. This corrected 
value of 1.001 ± 0.020 GeV is in good agreement 
with the theoretically corrected value from pion 
electroproduction of 1.014 ± 0.016 GeV. 

From PCAC, the pseudoscalar form factor Fp 
is predicted to be 



F P (q 2 ) = 



2M 2 F A (q 2 ) 
M%-q 2 ■ 



In the expression for the cross section, Fp(q 2 ) is 
multiplied by (m;/M) 2 . Therefore, in muon neu- 
trino interactions, this effect is very small except 
at very low energy, below 0.2 GeV. The effect is 
larger, about 5%, for tau neutrino interactions. 

3. UPDATED FORM FACTORS 

We have used an updated fit to the proton elec- 
tromagnetic form factors. The fit is similar to 
the one described in Ref. [5], but using a slightly 
different fitting function (described below), and 
including additional data to constraint the fit at 
low Q 2 values. Form factors can be determined 



from cross sections using the standard Rosenbluth 
separation technique |S] , which is sensitive to ra- 
diative corrections, or from polarization measure- 
ments using the newer polarization transfer tech- 
nique [S]. The polarization measurements do not 
directly measure the form factors, but measure 
the ratio Ge/Gm- Figures ^ El and [3] show the 
ratio of our fits divided by the dipole form, Gd- 
Figure 2| shows our fits to [i v G p E jG p M . The fit 
including only cross section data is roughly flat 
versus Q 2 (Q 2 — —q 2 ), while ratio decreases with 
Q 2 in the combined fit to cross section and polar- 
ization transfer data. Although the polarization 
transfer measurement is believed to have smaller 
systematic error, especially at high Q 2 , the ori- 
gin of this disagreement is not known. If this 
disagreement comes from radiative corrections to 
the electron, in particular two-photon exchange 
terms, then the polarization transfer extraction 
will give the correct ratio, but the overall scale 
of G P E at low Q 2 would be shifted down by «3%. 
Because the fit is constrained as Q 2 — » 0, there 
will not be an overall shift in G P E at low Q 2 , but 
there will be some uncertainty in the low Q 2 be- 
havior. Current experiments at JLab aim to bet- 
ter understand the source of the disagreement by 
looking at the recoil proton in elastic electron- 
proton scattering, thus minimizing the sensitivity 
to the dominant sources of uncertainty in previ- 
ous Rosenbluth separations. However, since this 
discrepancy is most prominent at high Q 2 , and 
the fit is constrained at low Q 2 , it has only a rel- 
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0.8548 
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-0.1287 


0.008912 
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CS 


3.226 


1.508 


-0.3773 


0.6109 


-0.1853 


0.01596 


CP 


CS 


3.188 


1.354 


0.1511 


-0.01135 


0.0005330 


-0.9005E-05 



Table 2 

The coefficients of the inverse polynomial fits for the G V B , G V M , and G^j. Fits using cross section data 
only, and using both cross section data and the Hall A polarization transfer data are shown separately. 
Note that these different polynomials replace Gd in the expression for G E , G P M , and G^. The first three 
rows of the table along with the fit of G P M Krutov et. al. 7 (see text) will be referred to as 'BBA-2003 
Form Factors'. 



atively small effect on the neutrino QE scattering 
cross section. 

To account for the fact that deviations from the 
dipole form are different for each of the different 
form factors, we fit electron scattering data for 
each of the form factors to an inverse polynomial 



Ge,m(Q 2 ) 



Gl M (Q 2 = o) 



1 + a 2 Q 2 + a 4 Q 4 + a 6 Q 6 + 

Table shows the parameters of our fit. We have 
done fits using cross section data only (for the 
proton) and fits using both both cross section 
data and polarization transfer data from JLab 
Hall A. For G^, the parameters in Table |21 are 
used for Q 2 < 6 GeV 2 . For Q 2 > 6 GeV 2 , the 
ratio of G V E /G P M is assumed to be constant: 

G E {Q)-G M {Q )^ (6 ^ 

Since the neutron has no charge, G E must be 
zero at q 2 —0, and previous neutrino experiments 
assumed G E (q 2 )=0 for all q 2 values. However, 
it is non-zero away from q 2 =0, and its slope at 
q 2 =Q is known precisely from neutron-electron 
scattering. At intermediate Q 2 , recent polariza- 
tion transfer data give precise values of G E (q 2 ). 
Our analysis uses the parameterization of Krutov 
et. al. 0: 



G E (Q 2 ) = -Mn 



1 + br 



G D (Q 2 ), 



4M 2 ' 



with a = 0.942 and b = 4.61. This parameteriza- 
tion is very similar to that of Galster et al. [5] , as 
shown in Figure 




Q 2 [GeV 2 ] 



Figure 4. Ratio of G P E to G P M as extracted by 
Rosenbluth measurements and from polarization 
measurements. The lines and symbols have the 
same meaning as Figure ^ 
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Figure 5 
is the Galster et al. fit [S] 
the Krutovet al. fit 0. 



The dashed line 
and the solid line is 
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Ratio, BBA 2003(CS+HallA,Krutov)/(Dip, GEn =0) 
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Figure 6. Ratio versus energy of predicted neu- 
trino (antineutrino) QE cross section using BBA- 
2003 Form Factors to the prediction using the 
dipole approximation with G E =0. 



The first three rows of Table [21 along with the 
fit of G E of Krutov et. al. [7j, will be referred 
to as 'BBA-2003 Form Factors'. For BBA-2003 
Form Factors, both the cross section and polar- 
ization data are used in the extraction of G E and 
ftp 

4. CROSS SECTIONS AND FITS TO M A 

Figure[S]shows the ratio versus neutrino energy 
of the predicted neutrino (antineutrino) QE cross 
section using our BBA-2003 Form Factors to the 
prediction using the Dipole Form Factors. The 
same comparison was performed using the pro- 
ton form factor extractions that included only the 
cross section data (i.e. excluding the polarization 
results). The results were nearly identical: the 
maximum difference between the cross sections is 
less then 0.3%. This is because the form factors 
extracted from the polarization transfer data and 
from the electron scattering cross section data are 
different only at high Q 2 , while the neutrino cross 
sections are mostly sensitive to the form factors 
at low Q 2 . 

There is a large difference in the predicted neu- 
trino (antineutrino) cross section between using 
the BBA-2003 Form Factors and using the Dipole 
Form Factor approximation. The difference is 3% 
at high energy and can become as much as 6% 
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Figure 7. Ratio versus energy of the predicted 
neutrino (antineutrino) QE cross section using 
from Krutov [x] to the prediction using 
0. In both cases, the dipole approximation 
for the other form factors. 
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Figure 8. Ratio versus energy of predicted neu- 
trino (antineutrino) QE cross section using our 
BBA-2003 Form Factor fits for G P E , G P M , and G% 
to the prediction using the dipole approximation. 
In both cases, we use G E =0. 



6 



at 1 GeV. As seen in Figure which shows the 



difference between GJ. 



Krutov and Gg = 0, 



all the difference at high energy and most of the 
the difference at low energy is due to G%. At 
the low energy region which is of interest for neu- 
trino oscillation experiments, both a non-zero Gg 
and the deviations from the dipole form are im- 
portant. This is also the case for the extraction 
of the axial form factor from neutrino data, since 
most of the neutrino differential cross section data 
are at low Q 2 . 

A 1% increase in either Ma or \qa\ increases 
the cross section about 1%. Replacing the old 
value of qa—— 1-23 with the more precise value of 
gA=— 1-267 increases the cross section by about 
2.5%. Using the more recent value of Ma of 1.02 
instead of the older value of 1.032 decreases the 
predicted cross section about 1%. F P has almost 
no effect on the cross section except at very low 
E v . Therefore, even a very conservative error 0] 
on Fp of 50% has very little effect. 

Previous neutrino measurements, mostly bub- 
ble chamber experiments, extracted Ma using the 
best known assumptions at the time. Changing 
these assumptions changes the extracted value of 
Ma- Hence, Ma needs to be updated using new 
form factors and up-to-date couplings. In this 
communication we will attempt to update the re- 
sults from three previous deuterium bubble cham- 
ber experiments. These are Baker et al. [H], Bar- 
ish et al. JU|> Miller et al. and Kitagaki et 
al. [T3|. Barish et al. and Miller et al. are the 
same experiment, with the analysis of Miller et al. 
including the full data set, roughly three times the 
statistics included in the original analysis. 

We start by calculating the shape of the Q 2 
distribution using the same form factors and cou- 
plings as used in the original extractions (includ- 
ing Ma). The flux is extracted from the flux 
figures shown in the original papers, which we 
parameterize using a spline fit. We extract the 
data and curves from their publications by pick- 
ing points off the plots, and fitting the points 
on the curves to a spline fit. These experiments 
did not use a pure dipole approximation, but in- 
cluded a correction to the dipole form as parame- 
terized by Olsson et al. [Tl|- They use gA=— 1.23, 



n -> p + p, , Baker 1981 
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Figure 9. Q 2 distribution from Baker et al. UJ. 
The dotted curve is their calculation taken from 
their Q 2 distribution histogram. The dashed 
curve is our calculation using their assumptions. 



+ n -> p + fj. , Barish 1977 
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Their calculation, M A =0.95 
Our calculation, M.=0.95 
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Figure 10. Q 2 distribution from Barish et al. |10j . 
The dotted curve is their calculation taken from 
their Q 2 distribution histogram. The dashed 
curve is our calculation using their assumptions. 



M y =0.84 GeV (yielding M2=0.7056 GeV 2 in- 
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p + (j. , Miller 1982 



Their calculation, M A =1.00 
Our calculation, M, = 1.00 
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Figure 11. Q 2 distribution from Miller et al. [TT] . 
The dotted curve is their calculation taken from 
their Q 2 distribution histogram. The dashed 
curve is our calculation using their assumptions. 



v + n -> p + fj, , Kitagaki 1983 



Their calculation, M, = 1.05 
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Figure 12. Q 2 distribution from Kitagaki et 
al. [T2|. The dotted curve is their calculation 
taken from their Q 2 distribution histogram. The 
dashed curve is our calculation using their as- 
sumptions. 



stead of M 2 =0.71 GeV 2 ), and a D 2 Pauli- 
suppression correction from Singh et al. 15 . Our 
calculations and their calculations are compared 
in Figures [SJ^H^TJ and 1 121 In these figures, our 
curves are normalized such that they match the 
previous curves at one point in the low Q 2 region, 
so that we can compare the Q 2 -dependence of the 
spectra. The y-axis is weighted events/bin (cor- 
rected for efficiencies). We reproduce the shape 
(Q 2 -dependence) of the calculations of Baker et 
al. liH, Barish et al. [TJ]|, and Kitagaki et al. |12| . 
but not Miller et al. ^JJ. As Miller et al. gives 
an updated result of Barish et al., they should 
be using the same calculation. Therefore, we do 
not understand the origin of the curve shown in 
Miller et al. 

Having reproduced the calculations under the 
same assumptions as in the original extractions, 
we perform our own extraction of Ma, using our 
calculation of the cross sections, but using the 
same input (gu and form factors) as assumed in 
the original extractions. Due to inefficiencies in 
reconstruction for very low Q 2 events, they do not 
use the first bin for fitting. Hence, we do not use 
the first bin for fitting or normalization either. 
We perform a binned maximum likelihood us- 
ing a formula from the Particle Data Group |16j . 
The experiments performed an unbinned likeli- 
hood fit, which we cannot reproduce since we do 
not have the individual events. Table EH gives the 
results of these fits. Using the same assumptions, 
we reproduce the fitted value of Ma from Baker 
et al., while disagreeing somewhat with the val- 
ues from Barish et al., Miller et al, and Kitagaki 
et al. The difference may come from the fact that 
we are forced to use a binned likelihood fit, rather 
than being able to reproduce their unbinned fit. 

Figures ITSI and ITU show the difference between 
using their value of Ma and our value of Ma for 
Kitagaki et al. and Barish et al. The plot for 
Kitagaki et al. appears to show that our value 
of M^^l.19 GeV is a better fit than their value 
of Ma- For Barish et al. its not clear which is a 
better fit. As previously stated, our Q 2 distribu- 
tion is slightly different than that of Miller et al. 
for the same Ma- Figure ITK1 shows their calcula- 
tion for their best fit Ma versus our calculation 
using our best fit Ma- The two shapes agree very 
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Our calculation, M. = 1.19 




Their calculation, M A =1.05 
Our calculation, M A =1.117 
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Figure 13. Q 2 distribution from Kitagaki et 
al. |12| . The dash curve is our calculation using 
our fit value of Ma=1.19 GeV. The solid curve is 
our calculation using their fit value of M^l.05 
GeV. 
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Our calculation, M, = 1.075 
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Figure 14. Q 2 distribution from Barish et al. |10|. 
The dash curve is our calculation using our fit 
value of Ma~1-Q75 GeV. The solid curve is our 
calculation using their fit value of Ma=1.01 GeV. 
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Figure 15. Q 2 distribution from Miller et al. 
The dotted curve is their calculation using their 
fit value of Af J 4=1.05 GeV. The dash curve is our 
calculation using our fit value of M^=1.117 GeV. 



well. Hence, we are able to reproduce the best fit 
shape of Miller et al. fit, but not with their value 
of M A . 

Finally, we extract Ma for each of these exper- 
iments using our calculations with the updated 
BBA-2003 Form Factors and gA value. Compar- 
ing this to our extraction with their input pa- 
rameters, we obtain change in Ma due to using 
updated values for gA and the form factors. Be- 
cause we sometimes obtain slightly different val- 
ues of Ma, even with the same input parame- 
ters, we take the difference between our extrac- 
tion with old and new form factors as the modifi- 
cation that should be applied to the previous ex- 
tractions, which were able to do a more detailed 
comparison to their data. Table gives the re- 
sults of these fits, which indicate that we should 
shift the value of Ma determined from deuterium 
by —0.025 GeV from the value quoted by these 
experiments. We also show that a shift in Ma of 
—0.050 GeV is required in going from the Dipole 
Form Factors to BBA-2003 Form Factors, keep- 
ing gA constant. If only cross section data were 
used for the form factor fits (Tabled rows 4 and 
5), the value of Ma would go up by 0.002 GeV, a 
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M A 
(published) 


updated Ma 
old params. 


updated Ma 
new params. 


AM A 
new-old 


AM a 
BBA-2003-Dipole 


Baker 1981 
Barish 1977 QJJ 
Miller 1982 [TTJ 
Kitagaki 1983 [11 


1.07 ± 0.06 
1.01 ± 0.09 
1.05 ± 0.05 


1.079 ± 0.056 
1.075 ± 0.10 
1.117 ± 0.055 

i.i94ig:iS 


1.055 ± 0.055 
1.049 ± 0.099 
1.090 ± 0.055 

1 - 1 '°-0.11 


-0.024 
-0.026 
-0.027 
-0.019 


-0.049 
-0.046 
-0.046 
-0.050 



Table 3 

Published and updated extractions of Ma (GeV) from deuterium experiments. The first value of Ma is 
the values extracted in the original publications. For Barish and Miller, we give their 'shape fit' value, 
since this value most closely reflects how we can calculate their Ma ■ The second value of Ma is from the 
analysis presented here, using the same input parameters (form factors and gA) as in the publications, 
while the third uses the updated parameters from tables |21 and ^ The last two columns show the change 
in Ma between the new and old input parameters, and the change when comparing the BBA-2003 and 
Dipole Form Factors (with gA fixed). 
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Figure 16. A comparison of the Q 2 distribution 
using 2 different sets of form factors. The data 
are from Baker et al. The dotted curve uses 
Dipole Form Factors with M A =l.lQ GeV. The 
dashed curve uses BBA-2003 Form Factors with 
M A =1.05 GeV. 



small effect. 

Figure [TBI shows the Q 2 distribution for Dipole 
Form Factors and Ma=1.10 GeV with the 
distribution for BBA-2003 Form Factors and 
M,4=1.050 GeV. When we modify the electro- 
magnetic form factors, the modification in Ma 
not only reproduces the original yield, it also re- 



produces the Q 2 distribution. Because there is no 
modification of the Q 2 dependence when strength 
is shifted between the electromagnetic and axial 
form factors, we conclude that the use of Dipole 
Form Factors will lead to an error in Ma of 0.050 
GeV, independent of the details of the experi- 
ment. 

Figures El an d El show the QE cross section 
for v and v using our most up to date assump- 
tions. The normalization uncertainty in the data 
is approximately 10%. We have used BBA-2003 
Form Factors, and have scaled down Ma from 
the old best fit of M A =1.026 ± 0.021 GeV to 
M J 4=1.00 GeV (which would have been obtained 
with the best vector form factors known today). 
The solid curve uses no nuclear correction, while 
the dotted curve EU uses a NUANCE gSJ cal- 
culation of a Smith and Moniz [23 based Fermi 
gas model for carbon. This nuclear model in- 
cludes Pauli blocking and Fermi motion, but not 
final state interactions. The Fermi gas model was 
run with a 25 MeV binding energy and 220 MeV 
Fermi momentum. The updated form factors im- 
prove the agreement with neutrino QE cross sec- 
tion data and give a reasonable description of 
the cross sections from deuterium. However, the 
data from heavy targets, including all of the anti- 
neutrino data, are systematically below the cal- 
culation, even with the Fermi gas nuclear cor- 
rections. The only experiment on heavy nuclei 
which agrees with the calculation is SKAT [T%| 
with neutrinos. The Fermi gas nuclear correc- 
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Figure 17. The QE neutrino cross section along with data from various experiments. The calculation 
uses Ma=1.00 GeV, £^=-1-267, M£=0.71 GeV 2 and BBA-2003 Form Factors. The solid curve uses 
no nuclear correction, while the dotted curve [21] uses a Fermi gas model for carbon with a 25 MeV 
binding energy and 220 Fermi momentum. The lower plot is identical to the upper plot with the E v axis 
limit changed to 2 GeV. The data shown are from FNAL 1983 12 , ANL 1977 [101, BNL 1981 0, ANL 
1973 [13, SKAT 1990 [Sj, GGM 1979 LSND 2002 [H|, Serpukov 1985 [HI, and GGM 1977 g2]. 
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Figure 18. The QE antineutrino cross section along with data from various experiments. The calculation 
uses Ma=1.00 GeV, g A =-l-267, M£=0.71 GeV 2 and BBA-2003 Form Factors. The solid curve uses no 
nuclear correction, while the dotted curve |2U uses a Fermi gas model for carbon with a 25 MeV binding 
energy and 220 MeV Fermi momentum. The data shown are from SKAT 1990 [Uj, GGM 1979 [25] . 
Serpukov 1985 21 , and GGM 1977 (22- 



tion may not be sufficient. Tsushima et al. |27j 
studied the effect of nuclear binding on the nu- 
cleon form factors. They stated that modifica- 
tions in bound nucleon form factors reduce the 
cross section relative to free nucleon form factors 
by 8%. We plan to study the nuclear corrections, 
adopting models which have been used in preci- 
sion electron scattering measurements from nu- 
clei at SLAC and JLab. In addition, we will be 
updating the extraction of Ma from other exper- 
iments, using the updated versions of the input 
parameters and electromagnetic form factors. 

Both the overall cross section and the in- 
dependence depend on the form factors chosen. 
However, even for input form factors that yield 
identical Q 2 -dependences, the cross section can 
differ significantly. Figure IT^l compares the cross 
section for two sets of form factors with identi- 
cal Q 2 distribution (shown in Fig. I16|) . The fig- 
ure shows that a cross section error as large as 
7.5% can occur from using this combination of 
incorrect form factors, even though they match 
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Figure 19. Ratio of cross section versus energy 
using BBA-2003 Form Factors with M A =1.00 
GeV versus Dipole Form Factors with M,4=1.05 
GeV. 
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Figure 20. Ratio of cross section versus energy 
using BBA-2003 Form Factors with Ma=1.00 
GeV versus Dipole Form Factors with Ma=1.11 
GeV. 



the l Q 2 shape' as well as the correct form factors 
do. The effects can be even larger than in this 
example, especially if the Q 2 -depcndence is not 
required to be identical. As shown by Y. Itow's 
at NUINT02 talk on the results from K2K 1 j, the 
initial K2K analysis (on a water target) assumed 
Dipole Form Factors. They obtained a value of 
Ma =1.11 GeV and saw a low Q 2 suppression with 
respect to the predicted distribution with Dipole 
Form Factors and Ma — 1.023. Figure I2TH shows 
the ratio vs. energy of cross sections predicted us- 
ing BBA-2003 Form Factors with M A =1.00 GeV 
to the prediction with Dipole Form Factors with 
Ma =1.11 GeV. The figure shows differences as 
large as 12%. This will clearly be important for 
future neutrino experiments at Fermilab and else- 
where: Even if the form factors which are used in 
the model are adjusted to describe the Q 2 dis- 
tribution measured in a near detector neutrino 
experiment, one would be predicting an incor- 
rect energy dependence of the QE cross section. 
Therefore, using the correct combination of form 
factors is important for determining the neutrino 
cross section and its energy dependence. An ac- 
curate measurement of QE scattering cross sec- 
tions requires an accurate measurement of both 
the normalized cross section versus energy as well 
as the shape of the Q 2 distribution. 



5. CONCLUSION 

We have made an updated extraction of nu- 
cleon electromagnetic form factors electron scat- 
tering, and have shown these have as much as 
a 6% effect for neutrino-nucleon QE scattering 
when compared with the standard dipole approx- 
imation. Inclusion of the new form factors yields 
to a reduction in the value of Ma extracted from 
neutrino scattering of 0.025 GeV. The agreement 
between the calculated neutrino and antineutrino 
free nucleon cross section and data is improved 
using the updated form factors, but is not spec- 
tacular, especially for data taken on heavy tar- 
gets. We have shown that for a fixed Q 2 depen- 
dence, the neutrino cross section and its energy 
dependence can be affected by more than 10% if 
the input shapes of the electromagnetic form fac- 
tors used are not correct. Hence, a complete un- 
derstanding of QE scattering scattering requires 
an accurate measurement of both the normalized 
cross section versus energy as well as the shape 
of the Q 2 distribution. In addition, nuclear ef- 
fects such as Pauli blocking and modification of 
nucleon form factors in bound nuclei need to be 
included. These effects are currently under in- 
vestigation. For example, the simple Fermi gas 
model can be modified to include a high momen- 
tum tail |2Hj or more sophisticated spectral func- 
tions. 

6. ACKNOWLEDGMENTS 

This work is supported in part by the U. S. 
Department of Energy, Nuclear Physics Division, 
under contract W-31-109-ENG-38 (Argonne) and 
High Energy Physics Division under grant DE- 
FG02-91ER40685 (Rochester). 

REFERENCES 

1. \protect\vrule widthOpt\protect\href -Qrttp : // 

2. Y. Fukada et al, Phys. Rev. Lett. 81 (1998) 
1562. 

3. C.H. Llewellyn Smith, Phys. Rep. 3C (1972). 

4. V. Bernard, L. Elouadrhiri, U.G. Meissner, 
J.Phys.G28 (2002), hep-ph[0107088]. 

5. J. Arrington, nucl-ex[0305009]. 

6. M. K. Jones et al, Phys. Rev. Lett, 84, (2000) 



1398 ; O. Gayou et al, Phys. Rev. Lett, 88 

(2002) 092301. 

7. A.F. Krutov, V.E. Troitsky, Eur. Phys. J. A16 

(2003) 285, hep-ph[0202183]. 

8. S. Glaster et al, Nucl. Phys. B32 (1971) 221. 

9. N.J. Baker et al, Phys. Rev. D23 (1981) 2499. 

10. S.J. Barish et al, Phys. Rev. D16 (1977) 3103. 

11. K.L. Miller et al, Phys. Rev. D26 (1982) 537. 

12. T. Kitagaki et al, Phys. Rev. D26 (1983) 436. 

13. T. Kitagaki et al, Phys. Rev. Lett. 49 (1982) 
98. 

14. M.G. Olsson, E.T. Osyporeski and E.H. Mon- 
say, Phys. Rev. D17 (1978) 2938. 

15. S.K. Singh, Nucl. Phys. B36 (1972) 419. 

16. Particle Data Group, Eur. Phys. J C15 (2000) 
196. 

17. W.A. Mann et al, Phys. Rev. Lett. 31 (1973) 
844. 

18. J. Brunncr et al, Z. Phys. C45 (1990) 551. 

19. M. Pohl et al, Lett. Nuovo Cimento 26 (1979) 
332. 

20. L.B. Auerbach et al, Phys. Rev. C66 (2002) 
015501. 

21. S.V. Belikov et al, Z. Phys. A320 (1985) 625. 

22. S. Bonctti et al, Nuovo Cimento 38 (1977) 
260. 

23. N. Armenise et al, Nucl. Phys. B152 (1979) 
365. 

24. G. Zeller, private communication. 

25. R.A. Smith and E.J. Moniz, Nucl. Phys. B43 
(1972) 605. 

26. D. Casper, Nucl. Phys. Proc. Suppl. 112 
(2002) 161. 

27. K. Tsushima, Hungchong Kim, K. Saito, 
nucl-th[0307013]. 

28. A. Bodek and J. L. Ritchie, Phys. Rev. D23 
(1981) 1070; ibid Phys. Rev. D24 (1981) 1400. 



